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Abstract 

In bipartite quantum systems of dimension 3x3 entangled states that are positive under partial trans- 
position (PPT) can be constructed with the use of unextendible product bases (UPB). As discussed in a 
previous publication all the lowest rank entangled PPT states of this system seem to be equivalent, under 
SL (g) SL transformations, to states that are constructed in this way. Here we consider a possible general- 
ization of the UPB constuction to low-rank entangled PPT states in higher dimensions. The idea is to give 
up the condition of orthogonality of the product vectors, while keeping the relation between the density 
matrix and the projection on the subspace defined by the UPB. We examine first this generalization for the 
3x3 system where numerical studies indicate that one-parameter families of such generalized states can be 
found. Similar numerical searches in higher dimensional systems show the presence of extremal PPT states 
of similar form. Based on these results we suggest that the UPB construction of the lowest rank entangled 
states in the 3x3 system can be generalized to higher dimensions, with the use of non-orthogonal UPBs. 



PACS numbers: 03.67.Mn, 02.40.Ft, 03.65.Ud 



INTRODUCTION 



One of the most striking features of quantum mechanics is found in the concept of entanglement. 
Considerable efforts have been made in understanding its properties and usefulness, with several 
different approaches taken [1]. One approach has been to study the geometrical structure of 
different convex sets of hermitian matrices that are related to subsets of quantum states with 
entanglement [2-7]. There are two convex subsets of the full set of density matrices, denoted V, 
that are particularly important in this discussion. One is the set S of non-entangled, or separable, 
states, and the other the set V of density matrices that remain positive under partial transposition 
with respect to one of the subsystems. The states of V are commonly referred to as PPT states. 
The non-entangled states are always PPT, so we have the set-theoretical relations S C V C V. 

In two previous papers we have examined the properties of entangled density matrices that 
are PPT, which means that they are contained in V but not in S [8, 9]. These states which are 
known to have bound entanglement are intrinsically interesting [10], but they are also interesting 
through the information they give about the set of entangled states. In fact, to establish whether a 
given density matrix is entangled or separable is usually a difficult problem, but only if the density 
matrix is PPT. This is so, since any state which is not PPT - a property that can easily be verified 
- is entangled. 

The states studied in [8, 9] were found by performing systematic numerical searches for PPT- 
states with specified ranks for the density matrix p and for its partial transpose p p . Such searches 
were performed for several bipartite systems of low dimensions, and in all cases, when the ranks 
were sufficiently low, the states were not only entangled but were typically extremal PPT states. 
An interesting result was that the lowest rank states of this type, with full local ranks, were found 
to have very similar properties, independent of their dimension. In particular they were all found 
to be extremal PPT states with a finite and complete set of product vectors in their kernel, and 
no product vectors in their image. 

For a bipartite system of dimension 3x3 there is a special subset of these lowest rank entangled 
PPT states that can be constructed by the use of unextendible product bases (UPB) [11]. The 
method provides a way to combine vectors in the two subsystems to form an orthogonal product 
basis that spans a five-dimensional subspace in the full Hilbert space, and which cannot be extended 
to include other product vectors orthogonal to the five states. The corresponding density matrix 
is constructed, up to a normalization, as the projection on the space orthogonal to the product 
vectors of the UPB. With the use of the numerical method described in [8] a large number of more 
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general entangled PPT states of rank 4 were generated, and it was shown that all of them, in a 
specific sense, were equivalent to density matrices that could be produced by the UPB construction. 
The equivalence classes were further shown to be parameterized by four real parameters, and the 
numerical results were interpreted as evidence for the conclusion that all entangled PPT states of 
rank 4 in the 3x3 system are covered by this parametrization. 

These results have motivated the present work, where we investigate a possible generalizion of 
the UPB construction which makes it applicable to systems of dimensions higher than 3x3. A 
simple copy of the UPB construction seems not possible, since in higher dimensions the orthog- 
onality requirement is too demanding. Instead we focus on other properties of the construction. 
In the 3x3 system the projection on the five-dimensional subspace spanned by the orthogonal 
product vectors of the UPB can be written as a sum of one-dimensional projections defined by the 
product vectors, and the corresponding density matrix is defined as the projection on the orthogo- 
nal subspace. The partial transpose of the density matrix takes the same form, when expressed in 
terms of a related, conjugate UPB. In the present paper we study density matrices of the similar 
projection form, but without the requirement of orthogonality. Such states can be found in the 
3x3 system, and we examine these states in some detail. We further focus on the question whether 
such generalized states can be found in systems of dimension higher than 3x3. 

The result is that we find such states in all the bipartite systems we have been able to examine. 
In addition to the 3x3 system this applies to systems of dimensions 3 x n with n taking values 
up to 6, and we have further examined the 4 x 4, 4 x 5 and 5x5 systems. The matrices we find 
are all of the projection form referred to above, with generalized, non-orthogonal UPBs in their 
kernel, and the partial transposed matrices are all found to have the same form. 

The organization of the paper is as follows. We first discuss in some detail the generalized UPB 
construction for the lowest rank PPT states of the 3x3 system. These states have rank 4 both for 
the density matrix p and its partial transpose p p , and we thus refer to them as (4,4) states. We 
illustrate the generalization by examining a special, symmetric case where the vectors of the UPB 
form a regular icosahedron, and find a one-parameter set of equivalent states that correspond to 
a linear deformation of the icosahedron. One of these states has an orthogonal UPB in its kernel, 
while for the general case the UPB is non-orthogonal. 

We next describe a numerical method to search for more general matrices, with less symmetric 
UPBs. The method specifies these density matrices p to be projections with rank equal to 4 and 
to have a positive partial transpose. All matrices that we find in this way are quite remarkably not 
only projections, but have the form we refer to in the generalized UPB construction. Furthermore, 
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the partial transpose p has the same form when expressed in terms of the associated, conjugate 
UPB. All density matrices found in this way are extremal PPT states. 

For the 3x3 system we have previously found numerically that generic extremal PPT states of 
rank (4, 4), found by the method described in [8], can be transformed by product transformations to 
a projection form with orthogonal UPBs in their kernels [9]. Here we further investigate whether 
such states can be transformed to a projection form with more general UPBs in their kernels, 
dropping the orthogonality requirement on the transformed UPBs. The results show that this is 
the case, and furthermore strongly indicate that the projections belong to one-parameter classes 
of equivalent projections, where density matrices defined by the orthogonal UPB construction are 
special cases. 

For the higher-dimensional systems the method we use to generate density matrices of projection 
form with specified ranks also works well. When the rank is chosen to coincide with the rank of 
the lowest rank extremal PPT states, as discussed in [8], we find matrices with the same properties 
as in the 3x3 system. The projection can thus be expressed in terms of a UPB in the kernel of 
the matrix, and the partial transpose has the same rank and structure when expressed in terms of 
its associated UPB. However, our further numerical studies indicate that a generic entangled PPT 
state of this rank cannot be transformed, by a product transformation, into a the projection form. 
In this respect the higher-dimensional systems seem to be different from the 3x3 system. 

We end with a summary and with a discussion of some of the questions that are left for further 
research. 

LOWEST RANK EXTREMAL PPT STATES IN THE 3x3 SYSTEM 

The construction of entangled PPT states of rank 4 in the 3x3 system, with the help of 
orthogonal UPBs, was originally discussed in [11]. A condition was there given for choosing a set 
of five vectors {4>k} in subsystem A and a corresponding set of vectors {xk} in subsystem B so 
that the set of product vectors {ipi- = <pk® Xk} would define an orthogonal UPB. From this set 
of product vectors a density operator of rank 4 could be constructed as a (normalized) projection 
operator in the following way 




(1) 



k=i 
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This form for p implies that the partial transpose p will have the same form when expressed in 
terms of a conjugate UPB, defined as {V>fc = 4>k ® X* k }i which is also a set of orthogonal product 
vectors. (The complex conjugation of the second factor implies that the partial transposition is 
performed with respect to subsystem B.) 

The important point is that the density matrix p, defined by the above expression, is necessarily 
entangled and PPT. It is entangled since there is no product state in the image of the matrix, and 
it is positive since it is proportional to a projection, which has only non-negative eigenvalues. For 
the same reason also the partial transpose is positive, and thus follows the PPT property of p. 
The density matrices defined by this construction form a proper subset of the entangled PPT 
states of rank 4 since the generic state of this type, according to our previous studies, will have a 
non-orthogonal UPB in its kernel [8]. 

The idea is now to relax the condition of orthogonality, but to keep other properties of the 
UPB construction. We examine this generalization first in the 3x3 system, but will subsequently 
examine the corresponding generalization for higher dimensional systems. The main condition is 
that the rank 4 density operator should be proportional to a projection operator, which we write 
as 

P=\{t-Q), Q 2 = Q (2) 

where Q is assumed to be of the form 

6 

Q = 5^2p k ip k ipl (3) 

with {V'fc = 4>k ® Xk} as a set of product vectors that defines a generalized UPB. The coefficients 
Pk define an unspecified set of real parameters, not necessarily all positive, with Ylk=iPk = 1- 

Note that the sum over the product vectors now runs from 1 to 6. The reason for this is the 
following. The five dimensional subspace spanned by the UPB will always include a total of 6 
product vectors [8]. When five of the product vectors are orthogonal, the 6th product state which 
is a linear combination of the other 5, is simply not included in the definition of p. We may view 
this as a special case of (2) and (3), with p$ = 0. However, when there is no subset of the product 
vectors that is orthogonal it seems natural to define the generalization so that all product vectors 
in the five dimensional subspace are included, as we have done above. 

The density matrix p defined by (2) and (3) is clearly entangled, since it has no product vector 
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in its image. Furthermore, the partial transpose p has the same form as p, when the product 
vectors ipk are replaced by the conjugate vectors ipk = 4>k ® X%i so that 

p p = \(l-Q p ) (4) 

with Q p given by 

6 

Q P = 5^p fc ^l (5) 
fc=i 

Also in this case the vectors ipk form a generalized UPB. However, it is not obvious from the above 
expressions that the rank of Q p is the same as the rank of Q, and neither that p p is positive. 
We only here note that in earlier numerical studies of the lowest rank extremal PPT states, we 
find that p and p p appear generically with equal ranks (in the 3x3 system as well as in higher 
dimensions) [8] . This implies that, for these states, if p is on projection form so is p p , as we shall 
show explicitly in a later section. In the following we shall simply assume that the generalized 
UPB construction will treat p and p p in a symmetric way with respect to the conditions (2) and 
(3), and thereby secure that p is both entangled and PPT. 

Since we lack an explicit prescription for choosing vectors <pk and Xk of the subsystems so that 
the product vectors ipk = <fik®Xk define density matrices p and p p that satisfy the above conditions, 
we have focussed instead on the question if we can, by numerical searches, find states that satisfy 
these criteria. This we have done first for the 3x3 system and subsequently for higher dimensional 
systems. In the 3x3 system we can, however, make an explicit construction of such states as a 
special case, and we shall discuss that case in the next section. 
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A SPECIAL CASE: THE ICOSAHEDRON 



An especially symmetric case, with non-orthogonal product states, is formed by the icosahedron 
construction described here. The vectors 0^ and Xk of the subsystems that combine into the product 
vectors of the generalized UPB are in this case all chosen to be real. The six vectors <pk define 
the six symmetry axes of a regular icosahedron that pass through its twelve corners. A particular 
choice of (non-normalized) vectors is specified by the following sets of Cartesian coordinates 

01 = (-0,0,1), 2 = (-1,0,0), 03 = (1,0,0), 

04 = (0,0,1), 05 = (0,-1,0), 6 = (0,1,0) (6) 

with = (y/E + l)/2 as the golden ratio. These are also the coordinates of six of the corners of 
the icosahedron, when the center is located at the origin. 

The second set of vectors, {xk}, is chosen as the same as set (6), but in order to form the correct 
combinations ipk = 4>k ® Xk: the 0^ and Xk vectors have to be differently ordered. A particular 
choice is Xk = 0(2fc+4)mod5 f° r k = 1,...,5 and %6 = — 06) but in total there are 60 acceptable 
orderings, related to the 60 rotational symmetries of the icosahedron. With all these orderings the 
six product vectors span a five-dimensional subspace of the Hilbert space. 

The vectors (6) define six equiangular lines, which means that the scalar product between all 
pairs of vectors are equal up to a sign. With normalized vectors the scalar products are 

^ = 0t0,=±-L, gB = X l Xl = ±-L, k + l (7) 

For the product vectors Vfc = 0fc ® Xk the scalar products are also equal up to a sign, and by 
choosing the particular ordering of the Xk, referred to above, we obtain 

9ki = 9ti9ki = -\ VMJ, M = l,...,6 (8) 
We now consider the following operator 

Q = lY,^4 ( 9 ) 
k=\ 

which is of the form (3) with all six coefficients equal, = 1/6. The condition that it should 
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define a projection, Q 2 = Q, can be written as 



= 5 ^2 9lk^l , 9ik = iptyk (10) 
and this should be satisfied for all k. When g^i is chosen as in (8) it gives the symmetric condition 

6 

5>* = ° (11) 

k=i 

and it is straight forward to check that this is equation is satisfied for the product vectors of the 
regular icosahedron. 

The product vectors {ipk} define a generalized UPB. Thus, there is no product vector orthogonal 
to the set, as can easily be checked, and the vectors are non-orthogonal. The set defines an entangled 
PPT state in the form of the density operator 

p=\^-\ii^i) (12) 

fc=l 

It is entangled since there is no product vector in its image, and it is PPT since the vectors are all 
real, and the partial transposition thus leaves the density operator invariant, p p = p. 
States that are related by non-singular product transformations 

(J = S P S^ , S = S A ® S B (13) 

where A and B refer to the two subsystems, all have the same characteristics as being separable or 
entangled. They also share the property of being positive or not under partial transposition and 
they have the same rank. This is trivially the case for unitary product transformations, but it is 
also true for non-unitary transformations, in which case the transformed matrix p' should further 
be normalized to unity. The operators Sa and Sb can be restricted, without loss of generality, 
to be unimodular, and for this reason we refer to the relation as SL (g) SX-equivalence (or simply 
SL-equivalence) [9]. For the density operator (12) there is a subset of such SX-equivalent states 
that can be written in the form of projections (3). Clearly unitary product transformations will 
leave this form invariant, but it is of interest to note that there is a one-parameter set of non- 
unitary transformations that also leaves the projection form invariant. As opposed to the unitary 
transformations these transformations will change the coefficients in the expansion (3). 
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FIG. 1: The icosahedron construction. A regular icosahedron, shown to the left, defines a set of six equian- 
gular lines through the twelve corners. A corresponding set of six vectors along these lines are shown in the 
figure. Tensor products of six pairs of such vectors define the vectors of the UPB, as explained in the text, 
and all of these appear with equal weight in the construction of the corresponding entangled PPT state. To 
the right a stretched icosahedron is shown where five of the vectors are used to define an orthogonal UPB. 
These vectors point along edges of a regular pyramid with a pentagon as base, which appears upside down 
in the figure. This construction has previously been referred to as the Pyramid [11]. The corresponding 
UPB defines an entangled PPT state where the five product states appear with equal weight, and where the 
sixth vector does not appear. 

The non- unitary S'L-transformations are of the symmetric form S = Sa &) Sb with Sa = Sb, 
where Sa rescales the length in the direction of one of the symmetry axes of the icosahedron. In the 
following we choose this to be the direction of 06, so that the transformation squeezes or elongates 
the icosahedron along this axis. This can be described in terms of a streching parameter A, so that 
the unit vectors pointing towards the corners of the deformed icosahedron are 

0fe(A) = N(X) + (A - l)(4>l<f>k) 06 ) j M6; 6 (A) = 06 (14) 

with N(X) as a normalization constant. For symmetry reasons, the coefficients Pk(X) are all equal 
for k 6. It is straight forward to show that 

4 + 2A 2 -A 4 (1~P 6 (A)) 
Pe{X) = 20 + 10A2 ; PkW = 5 V M 6 (15) 

The stretching parameter is chosen with A = 1 for the regular icosahedron with p^ = 1/6 
for k = 1,...,6, in which case all product vectors contribute equally. When A = the squeezed 
icosahedron collapses to a plane, at which point the set of product vectors cease to be a UPB 
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since the vectors (f>k for k = 1 , . . . , 5 all lie in the two-dimensional plane orthogonal to <pQ . For the 
particular value A = y/2cj) the five product vectors, with k = 1,2,. ..,5, become orthogonal, and 
thereby define an orthogonal UPB. This is the Pyramid construction discussed as a particular case 
of a UPB construction in [11]. In this case the coefficients are pu = 1/5 for k / 6 and p$ = 0. The 
relation between this set and the stretched icosahedron is illustrated in Fig.l. 

THE GENERAL CASE OF RANK (4, 4) ENTANGLED PPT STATES 

We focus now on general entangled PPT states of rank (4, 4) in the 3x3 system, which are 
necessarily also extremal PPT states. In a previous study [8] we have numerically produced a 
large number of such states, and we have found that they all can be transformed, by product 
transformations, to states of the special form given by (1), with an orthogonal UPB in the kernel 
of the density matrix [9]. The question that we now discuss is whether these are only special cases, 
and that more general states on the projection form given by (2) and (3) can be found. In order 
to examine this we have applied two different numerical methods to search for density matrices of 
the desired form. 

The first method applies search criteria which do not directly refer to the conditions (2) and 
(3). Instead the search is for density matrices p with the correct rank 4, which are PPT and of 
projection form. This means that all the non-vanishing eigenvalues should be equal. The method is 
essentially the same as used in [8] to search for PPT states of rank (4, 4). Here however, we impose 
no restriction on the rank of p p , only that det p p > 0. It is a linearized, iterative method which 
specifies a certain number of the eigenvalues of p to vanish, and in this case also the remaining 
eigenvalues to be equal. We refer to [8] for details concerning the iterative method. A similar 
approach is also used in the second method, to be discussed below. 

The result is that the numerical method works well, and a large number of states which satisfy 
the search criteria have been found. Even if we do not introduce any explicit constraint on p p , 
only the conditions that p is PPT and is proportional to a projection, we find always p p also to 
be proportional to a projection, so that both p and p p satisfy the condition (2). Futhermore we 
find in all cases that p (and p p ) satisfy condition (3), when expressed in terms of the vectors of 
the UPB, even if that is not included as a condition in the searches. 

To gain some more information about the type of density matrices we want to obtain, we have 
applied a second numerical method. It introduces a search for a product transformation which 
maps, if possible, a chosen rank (4, 4) extremal PPT state into the form of a projection. We do 
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not, in this search either, impose any condition on how this projection is expressed in terms of 
product vectors. Also this method applies a linearized, iterative approach to determine the product 
transformation, which we now outline. 

Let us denote the initial density matrix by po and the final density matrix by p. Their relation 
should then be of the form 

p = S Po S^ , S = S A ® S B (16) 

and the condition that the transformed density matrix is a projection gives 

SpoS^SpoSi = S Po S^ (17) 

(Note that the density matrix p is then not normalized to unity.) The transformation S is assumed 
to be non-singular, and the condition can then be simplified to 

Po&Spo = p (18) 

The unknown transformation matrix S we parametrize by expressing Sa and Sb as linear 
combinations of a complete set of 3 x 3 hermitian matrices. The corresponding product transfor- 
mations are parametrized by 17 real parameters p^, which we interpret as components of a vector 
/Li. Furthermore we express the rank 4 matrices as 16-component vectors in the space of hermitian 
operators within the image of po- Written as a vector equation, Eq. (18) has the form 

F(fjL) = po (19) 

Assume now that fx' is a trial vector that gives an approximate solution to Eq. (19). We write 
the deviation from the true solution as A fx = fx — fx', and treat this as a perturbation. To first 
order in Afx the equation reads in matrix form 

BAfx = p - F (20) 

with B as a real, non-quadratic 16 x 17 matrix with elements Bk n = dFk/dp n , and where both 
F and B are evaluated with the trial vector fx' . By multiplying with the transposed matrix B T 
and introducing the positive, real symmetric matrix A = B T B, as well as a = B T (p — F), the 
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equation can be written in the form 

AA/x = a (21) 

where a and A are determined by the trial vector /x' and A/i is the unknown to be determined by 
the equation. Written in the form (21) the equation is well suited to be solved numerically by the 
conjugate gradient method [12]. 

An iterative approach is used to find a solution of the original problem. A starting point 
H 1 is chosen for the trial vector, and this determines the initial versions of a and A. Eq. (21) 
is then solved numerically to give a first solution A/^. The trial vector is then updated with 
fi 2 = Mi + A/^, this vector is used to improve a and A, and a new improved solution A/x 2 of (21) 
is found. If repeated iterations of this procedure leads to convergence, in the sense A/z — > 0, the 
limit value of /z gives a solution to the original problem (19) and (18). 

A possible problem with this method is that the solution we find may correspond to a singular 
transformation matrix S, which will not give a density matrix p with correct properties. However, 
for the 3x3 system, by repeatedly applying the method with different starting points, we have 
found that usually the convergence of the method is quite rapid, and the solution that we find 
corresponds to a non-singular product transformation. 

The result from applying the above method to the rank (4, 4) states of the 3x3 system is that 
for a large number of initial density matrices that we have used, we can in all cases transform the 
density matrix by a non-singular product transformation to a matrix with the form of a projection. 
In all cases the density matrices p that we find have the form specified by (2) and (3), when 
expressed in terms of the (non-orthogonal) UPB in the kernel of the density matrix. Furthermore, 
also the partially transposed density matrices p p are, in all cases, proportional to projections, and 
can be written in the form given by (2) and (3). This happens even if only p is required to have 
the form of a projection in the numerical search. 

By keeping po in (18) fixed and choosing repeatedly different starting points fi 1 in (19) for 
the first iteration, we get different results for the the density matrix p determined by the above 
method. This indicates that there is a large number of different states on projection form within 
a given set of SL-equivalent (4, 4) states. In fact there are several indications that there is a one 
parameter set of such states, which can be related by non-unitary product transformations, within 
any given equivalence class of extremal PPT states with this rank. The first indication is simply 
based on parameter counting. As already discussed the number of equations that determine a 
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product transformation that transforms a state of rank 4 to the projection form is 16 while the 
number of parameters to be determined by the equations is 17. 

The second indication is found by listing the parameter sets {pk, k = 1, 2, 6} for 5L-equivalent 
states on projection form that we generate by our method. As shown in Fig.2a a large set of different 
distribution is found, consistent with the assumption that there is a continuum of equivalent states. 
If one of the parameters is restricted to a very limited interval, the corresponding sets {pk} seem 
either to be identical (up to the limitation set by the fixed parameter), or to divide into a small 
number of distinct groups, each of which are essentially identical. This is illustrated in Fig. 2b. This 
is consistent with the assumption that the full set is specified, up to a discrete set of possibilities, 
by a single continuous parameter. 

The third indication is that we have checked numerically, in a few examples, that if po in (18) is 
already on projection form, then this equation has solutions of the form S = (l^+eL^g^ls+eLe) 
with e an infinitesimal parameter. We find in each case that this set of infinitesimal product 
transformations is uniquely determined if we require S to be hermitian and not unitary. 

All this seems to show that the picture is similar to that of the icosahedron case, so that any 
extremal PPT state of rank (4, 4) is SX-equivalent to a one parameter set of states on the projection 
form given by (2) and (3). And for each of these the partial transposed density matrices has the 
same projection form. (Note that the trivial equivalence under unitary product transformations is 
not included in this parameter counting.) 

HIGHER DIMENSIONS 

The lowest rank extremal PPT states in the 3x3 system have their counterparts in higher 
dimensional systems. Thus, in previous numerical studies of PPT states in several systems of 
dimension Na x Nb, with Na and Nb larger than 2, we have found such states with rank given by 

r = N A + N B -2 (22) 

for both p and p p [8]. These PPT states are typically extremal, all with a UPB with a finite 
number of (non-orthogonal) product vectors in their kernel. This is a situation which is different 
from what we find for entangled PPT states with other values of the rank. Based on these results 
we have conjectured that, quite generally for higher dimensional systems, the number (22) gives 
the lowest rank of extremal PPT states with full local ranks for a bipartite system of dimension 
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FIG. 2: a) The figure shows the distribution of the smallest coefficient in the set {pk} for approximately 9500 
SX-equivalent states on projection form, all of which are SX-equivalent to a generic rank (4, 4) extremal 
state. We observe that we find a large distribution that seems to cover all possible values from up to 
an upper limit which is smaller than 1/6, which is the upper limit of the icosahedron states. This implies 
that such generic rank (4, 4) extremal states are less symmetric than the icosahedron-class of states. The 
continuous nature of the distribution also supports the one-parameter family theory, while we do not have 
an explanation for the particular profile of the distribution, b) This plot shows some of the sets {pk} where 
the smallest parameter is specified within three very limited intervals. With this limitation the rest of the set 
is found to be uniquely defined within essentially the same limitation. In one case, however (the dashed blue 
curve), there are two distinct groups of parameter sets with the same smallest coefficient. These obervations 
support the claim that the SL-equivalent states on projection form defines a continuous one-parameter set 
of states. 

Na x Nb- The number of product states for a UPB associated with any of these states is [8] 



(N A + N B -2)\ 
(N A -l)\(N B -l)\ 



(23) 



and the number of linearly independent product states is 



d = N A N B -N A -N B + 2 



(24) 



In Table 1 we have listed the relevant number of product states for the systems which we have 
studied numerically in this work. The table also includes a list of dimensions for the sets of PPT 
states with the relevant ranks in these systems, as well as the number of parameters needed to 
specify the classes of SX-equivalent matrices of this type. These numbers are based on numerical 
studies that we have now performed, and which are described below. 

The method used to determine the dimensions of the sets of PPT states is based on the counting 
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system 


ranks 


p/d 


dimensions 


3x3 


(4.4) 


6/5 


36/4 


3x4 


(5,5) 


10/7 


55/9 


3x5 


(6,6) 


15/9 


77/13 


3x6 


(7,7) 


21/11 


102/16 


4x4 


(6,6) 


20/10 


75/15 


4x5 


(7,7) 


35/13 


96/18 


5x5 


(8,8) 


70/17 


119/23 



TABLE I: Ranks and numbers of product states in the kernel of lowest rank extremal PPT states (with full 
local rank) in a series of bipartite systems of low dimensions. The first column gives the dimensions of the 
two subystems, the second column the ranks of the density matrices and their partial transpose (which are 
equal for all these states), the third column gives the total number of product states in the kernel of the 
density matrices as well as the number of linearly independent product vectors, and the fourth column gives 
the dimension of the set of the density matrices and the number of parameters needed to parameterize the 
classes of SX-equivalent density matrices. 

of different ways to make small perturbations away from a density matrix p of a given rank (m, n), 
in such a way that the ranks of the matrix and its partial transpose are preserved. Consider then 
po to be an extremal PPT state, and (m, m) to be the (symmetric) ranks of the density matrix 
and its transpose. A perturbation of this state we write as 

p = p + ea, e << 1 (25) 

Further, let P be the projector onto the image of po and Q the projector onto the image of p p . 
The following conditions, (/ — P)a(I — P) = and (/ — Q)a p (I — Q) = 0, secures that the ranks 
(m,m), to first order in e, do not change under the perturbation. We rewrite the conditions as 

a = Pa + aP - PcrP = <j> P {a), a p = Qa p + a p Q - Qa p Q = 4>q((t p ) (26) 

The maps <pp and 4>q define linear operators P and Q that act as projectors on the real vector 
space of hermitian matrices. The previous equations can therefore be written as 

Pa = a, Q(T P = a p (27) 

with a and <j p viewed as vectors. The partial transpose can further be expressed as a linear 
operator II in this space, so that II <r = a p [13]. With the definition Q = IIQII, the two 
equations in (27) can be combined in the single equation 

PQPa = a- (28) 



15 



and the number of linearly independent matrices a that obey this equation will then be precisely 
the dimension of the set of matrices with rank (m, m) in which po sits. 

The number of linearly independent matrices a can be found by calculating the matrix PQP 
and counting the number of its eigenvalues that are equal to 1. As we are interested in the 
dimension of the set of extremal density matrices, with the specified rank, we have also checked 
whether the states close to po are typically extremal. This has been done by repeatedly perturbing 
po in different directions and checking the perturbed states for extremality. Since the rank, for a 
finite perturbation, will generally increase as a higher order effect in e, we have in each perturbation 
corrected for this by including a second search, if necessary, for a closeby density matrix with correct 
rank, before checking for extremality. 

The result is that for all the states listed in Table 1, we find that the density matrices with the 
correct ranks (m, m) in the neighborhood of a chosen extremal density matrix are typically also 
extremal. This we take as a clear indication that the dimension we calculate is also the dimension 
of the set of extremal states with the given rank. 

To determine from this dimension the number of parameters that is needed to parametrize 
the classes of SX-equivalent states, we have subtracted the number of parameters that specify a 
product transformation of the form SL(Na, C) x SL(Nb, C). Each factor is specified by 2N 2 — 2 
real parameters, and the total number of parameters is therefore 2(iV| + Ng) — 4. The numbers 
given in Table 1 are obtained by such subtractions. In particular the number of parameters found 
in this way for the (4, 4) states of the 3x3 system is 4. This agrees with the conclusion reached 
in [9]. 

PROJECTION OPERATORS IN HIGHER DIMENSIONS 

The same numerical methods that have been used to study the (4, 4) states of the 3x3 system 
we have also applied to the extremal PPT states in higher dimensional systems. A specific case is 
the 4x4 system, where the relevant states are of rank (6, 6). These density matrices have a UPB 
with 20 product vectors in their kernel, 10 of these being linearly independent. As shown in Table 
1 we find that these states form a 75-parameter subset within the 255-parameter set of density 
matrices of the 4x4 system. The number of parameters determining a product transformation 
in the 4x4 system is 60, which leaves us with 15 parameters with which to parameterize the 
equivalence classes. 

The method we use to search for PPT states of the correct rank and on projection form, works 
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well also for this system, and we have by use of the method generated a large number of density 
matrices with these characteristics. We find, precisely as in the 3x3 system, that the density 
matrices always have a diagonal form, similar to (3), when expressed in terms of the product 
vectors of the UPB. In the 4x4 system the correct expressions for the density matrix p and the 
corresponding projection Q is 

1 20 
p=-(l-Q), Q = 10 J^PkMl, Q 2 = Q (29) 

k=i 

with tpk = 4>k <8> Xk as the product states of the UPB. Precisely as in the 3x3 system we also here 
find that the partially transposed operator Q p is a projection with the same rank as Q, so for the 
density matrices that we find there is complete symmetry between p and p p , 

1 20 

p p = -(t-Q p ), Q p = ioJ> fc ^I, (Q P ? = Q P (30) 

k=i 

where 4>k = <Pk® Xk are the product vectors of the conjugate UPB. 

The states we find have generally pk / for all k = 1,2, ...,20. However, while most of these 
coefficients are positive, a small number of them will usually be negative. This is different from 
what we find in the 3x3 system, where the typical situation is that all pk are positive, but where 
we occasionally find one of them to be negative. 

Also for the 4x4 system we have investigated the possibility of transforming a generic extremal 
PPT state of rank (6, 6) to the projection form. The method we use is the same as for the 3x3 
system, where we search for solutions to Eq. (18), with po as a density matrix with the correct 
rank, which is generated by the method described in [8]. However, as opposed to the case with 
the (4, 4) states of the 3x3 system, we have only for special choices of matrices po been able to 
find product transformations that transform po to the projection form. In the general case the 
iterative method that we use does not converge to an acceptable solution. Instead it shows a slow 
convergence towards a singular transformation matrix. 

The lack of convergence in this case may be a consequence of the increase in the number of 
variables in the problem, and therefore to a a decrease in the efficiency of the iterative procedure. 
But a clear possibility is that in the 4x4 system a generic extremal PPT state of rank (6, 6) 
cannot be transformed by a product transformation to the form of a projection. In fact, the form 
of the equation (18) that we seek solutions for may indicate that this is the case. By counting the 
variables of the equation we find that the set of equations is underdetermined in the 3x3 system, 
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but it is overdetermined in the 4x4 system as well as in other higher dimensional systems. (That 
does not, however, exclude the possibility that for the special density matrices that we consider 
there should exist solutions to the equation.) 

The other higher-dimensional systems that are listed in Table 1 have been studied by the same 
methods as the 3x3 and 4x4 systems, and the results are essentially the same as for the 4x4 
system. This means that the searches for PPT states p with rank m as specified in Table 1, which 
are proportional to projection operators, are in most cases successful. By varying the initial value 
in the search we have therefore been able, for most of the listed systems, to generate a large number 
of different solutions. For systems of dimension N = NaNb > 20 the iterative methods become 
rather slow, so the number of solutions we have found for them is somewhat smaller, though we 
find solutions also there. In all cases we find density matrices of the same form as shown in (29). 
We also find in all cases that the partially transposed matrix p p is of the same rank as p. It is also 
a projection, and therefore we have a complete symmetry between p and p p , similar the one given 
by (29) and (30) for the 4 x 4 system. 

For comparison we have also made searches for PPT states that are proportional to projections 
with ranks higher than those indicated in Table 1, in which case there is typically no UPB in the 
kernel of the matrix. The result is that we are able to find also such density matrices, but now the 
situation is different. In these cases the typical solution p has a partial transpose p p with a higher 
rank and which is not proportional to a projection. 

We have for all the listed higher dimensional systems also performed searches for product 
transformations that transform a generic extremal PPT state, with the specified rank, to projection 
form. For these systems the result is the same as for the 4x4 system, that the searches in most 
cases are unsuccessful. Thus, only for the 3x3 system do we find that we are able to transform 
the generic extremal PPT states of the given rank into the projection form. 

We summarize the main findings for the higher dimensional systems in the following way. For 
all the bipartite systems we have studied we have been able to generate entangled PPT states on 
projection form, with both p and p p satisfying the conditions given by (2) and (3). Whereas in the 
3x3 system all the entangled PPT states with the relevant rank seem to be SL equivalent to the 
states on projection form, in higher dimensions these states seem instead to form a proper subset 
of the entangled PPT states with the given rank. In addition to these results we have made other 
observations that may be relevant for the generalized UPB construction, and we discuss some of 
them in the next section. 
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STATES ON PROJECTION FORM AND REAL INVARIANTS 



We first return to a point briefly mentioned earlier in the paper about a difference in the UPB 
construction which appears when orthogonal UPBs are replaced by generalized UPBs. In the 
former case the PPT property follows directly from the projection form of p, since this implies also 
p p to be proportional to a projection, and therefore to be positive. With the use of a generalized 
UPB, we have no proof for this to be the case. Clearly the UPB construction will be simpler if we 
only have to take into account conditions that apply to p and not to both p and p p , and for that 
reason we shall examine this point a bit further. 

We first demonstrate that symmetry between p and p p in the above respect is directly linked 
to a condition of equal ranks for the two matrices. At this point we make no specific assumptions 
about the matrix p except for the normalization Tr[p] = 1. We refer to the rank of the matrix as 
r and P as a projection on the image of p. We have 

Tr[(p Tr[p 2 }- \ (31) 

which implies 

Tr[p 2 ] > I (32) 

with equality when p = P/r. We next assume precisely this to be the case, so that p is proportional 
to a projection and therefore Trip 2 ] = V r - We further introduce s as the rank of p p and apply 
the same argument to this matrix, so that 

Trl(p p ) 2 ] > J (33) 

On the other hand, the trace of the matrix squared is invariant under the operation of partial 
transposition, Tr[(p p ) 2 ] = Tr[p 2 ] = 1/r, and this implies the inequality 

s > r (34) 

This means that the rank of p p is larger or equal to the rank of p when the latter is proportional 
to a projection, and there is equality between the ranks if and only if p p is also proportional to a 
projection. From this we conclude that p p satisfies the same projection condition (2) as p if and 
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only if they have equal ranks. 

As a particular application let us assume that p is of the projection form given by (2) and 
(3), and that it can be transformed to a situation where the product vectors ipk = 4>k® Xk of the 
generalized UPB are decomposed into vectors <pk and Xk of the subsystems where the vectors of 
one of these sets are real. This may be taken as subsystem B, so that x% = Xk- As a consequence 
the density matrix is invariant under partial transposition, and therefore the ranks of p and p p are 
equal, not only after the transformation, but also in the reference frame where p is on projection 
form. This is so since the rank is invariant under the non-singular product transformations. As 
follows from the discussion above the partial transpose p p will in the same reference frame also be 
proportional to a projection, and hence positive. 

The condition of real vectors is an obvious way to secure the density matrix and its partial 
transpose to have the same rank. On the other hand, this condition could possibly be too restrictive 
to make it interesting. However, in [9] a set of SL&SL invariants was introduced that characterize 
the product vectors of the generalized UPBs, and for the relevant rank (4, 4) states of the 3x3 
system they were found to be real. This is a non-trivial result, since for a UPB more generally 
the invariants will be complex. When all the invariants are real it follows that the vectors of the 
UPB can be transformed to real form. For the higher dimensional systems the situation is more 
complex. The 3 x n systems (n > 3) that are available for numerical study we have found all to 
have real invariants for the set of vectors associated with one of the subsystems, but complex for 
the other subsystem. However, for other systems, like the 4x4 system, the generic lowest rank 
extremal PPT states are characterized by SL invariants that are complex for both subsystems. 

A similar study of the states referred to in the present work, where the density matrices are 
on projection form, however shows a somewhat different picture. For all the systems we have 
examined we find the invariants of the states that have been generated are all real for at least 
one of the subsets of vectors {4>k} or {xk}- Thus, in all these cases one of the sets of vectors 
can be transformed to real form. This indicates a close connection between the reality condition 
referred to above and the condition of projection form for p and p p . Note however, that to impose 
the reality condition and the projection condition simultaneously may be too restrictive, since the 
conclusion we can draw from our studies is only that the density matrices we have found, which 
satisfy the projection condition, can be transformed to a form where one of the subsets {4>k} or 
{Xk} consists only of real vectors. 
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CONCLUSIONS 



The motivation for the present work has been to examine the possibility of generalizing the 
method of constructing entangled PPT states with the use of unextendible product bases (UPB). 
The established form of this construction is to use a set of orthogonal product states, with no 
product states in the orthogonal subspace, and to define the corresponding density matrix as a 
projection operator with the UPB in its kernel. The method applies particularly to rank (4, 4) states 
in bipartite quantum systems of dimension 3x3, and as previously shown in numerical studies all 
entangled PPT states with this rank seem to be equivalent under product transformations to the 
states constructed in this way [9] . 

In higher dimensional bipartite systems there are states that share many of the properties with 
the (4, 4) states of the 3x3 system. They have all a finite and complete set of product vectors 
in their kernel and no product vector in their image, and they all seem also to share the property 
of being the lowest rank extremal PPT states in the system under consideration. The idea is to 
extend the UPB construction to these states. It seems not to be possible to extend the construction 
with orthogonal product vectors to higher dimensions, and we have focussed on the possibility of 
using non-orthogonal UPBs which keep some of the other properties of the original construction. 
Our assumption is that we can define these density matrices more generally as projections which 
can be expressed in a particular way in terms of product vectors of the non-orthogonal UPB. 

We have first examined this generalization for rank (4, 4) extremal PPT states of the 3x3 
system. In this case an explicit example can be found, with the vectors of the generalized UPB 
defined by the symmetry axes of an icosahedron. All the six product vectors of this UPB appear 
with equal weight in the definition of the corresponding extremal PPT state. We can further show, 
by a linear deformation of the icosahedron, that a one-parameter set of different density matrices 
on projection form exists, where the matrices are equivalent under product transformations. 

To study more general states we have applied numerical methods. The first method is to search 
for PPT states which have rank 4 and which are proportional to projections. The searches have 
been used to generate a large set of such states and these states have been found always to be on 
the form suggested by the generalized UPB construction, where p and p p have equal ranks, are 
both proportional to projections and have the same form when expressed in terms of the product 
vectors of the UPBs associated with the two matrices. 

A second method has been used to check whether generic extremal PPT states of rank (4, 4) 
are equivalent under product transformations to density matrices of the form suggested for the 
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generalized UPB construction. This has been demonstrated numerically and the numerical studies 
suggest that each state is equivalent to a one-parameter set of density matrices of this form where 
the matrices associated with an orthogonal UPB constitute a discrete subset. 

To examine the relevance of the generalized UPB construction in bipartite systems of higher 
dimensions, we have first studied numerically the dimensions of the sets of extremal PPT states 
of the corresponding ranks. We have then applied the same numerical methods as used for the 
3x3 system to search for density matrices of the suggested form. The search is thus for PPT 
states on projection form, with the specified rank, and the result are similar to those found in the 
3x3 system. For all the systems a large number of states which satisfy the search criteria have 
been found, and they all show a complete symmetry between the density matrix and its partial 
transpose. They have the same rank, are both projections and can be expressed in terms of the 
product vectors of the associated UPBs in the same way. This demonstrates that the requirements 
suggested for a generalized UPB construction are satisfied for a large set of states in all these 
systems. 

We have however observed a difference between the 3x3 system and higher-dimensional systems 
when performing searches for states of projection form that are equivalent under product transfor- 
mations to a randomly generated extremal PPT state of the correct rank. As we mentioned, in the 
3x3 system the results indicate that any extremal PPT state with the specified rank is equivalent 
to a one-parameter set of states on projection form, where the states associated with orthogonal 
UPBs form a discrete subset. In higher dimensions similar searches have been unsuccessful in 
most cases. This may suggest that the density matrices in higher dimensions that are of the form 
specified in the searches form a proper subset of the full set of extremal PPT states with the given 
rank. This is further supported by the fact that the states on projection form of the correct rank 
in higher dimensions are found to have real SL- invariants, in contrast to the all complex invariants 
of randomly generated states of the same rank. 

Based on these results we suggest that the generalized UPB construction may be relevant for 
construction of low-rank extremal PPT states in higher dimensions. Our study thus shows that 
extremal PPT states of the suggested form exist in all the higher-dimensional systems we have been 
able to examine. However, we still lack a concrete prescription for constructing generalized UPBs 
which define density matrices with the right properties. This problem, to find such a prescription, 
and also the question about how general the density matrices of this form are, we shall therefore 
have to leave as interesting questions for further studies. 
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